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Abstract 


This  study  aims  to  formulate  a  method  of  predicting  the 
far  field  pattern  of  a  parabolic  dish  antenna  with  two 
moveable  flat  plates  mounted  symmetrically  on  either  side 
of  the  feed  horn.  The  approach  taken  has  been  to  first 
analyze  the  radiation  pattern  of  the  antenna  with  the  disks 
at  certain  heights  out  from  the  surface  of  the  dish.  To  do 
this  the  near-field  radiation  in  amplitude  and  phase  was 
measured  over  a  plane  surface  in  the  near-field  and  the 
values  were  then  transformed  into  the  far  field  using  a 
Fast  Fourier  Transform. 

Far  field  pattern  values  of  the  antenna  were  directly 
measured  for  each  setting  of  the  plates.  The  results 
obtained  from  the  Fast  Fourier  Transform  of  the  near  field 
data  were  in  good  agreement  with  the  values  obtained  by 
measurement . 

Finally,  an  approximate  model  of  the  antenna  was 
developed  and  implemented  as  a  computer  program.  This 
model,  while  relatively  unsophisticated,  provided  some 
insights  into  the  changes  in  the  near  field  phase 
distribution  caused  by  the  moveable  circular  flat  plates. 


MODIFICATION  OF  PARABOLIC  DISH 


ANTENNA  PATTERN  USING  SYMMETRICALLY 
PLACED  CIRCULAR  FLAT  PLATES 


I .  Introduction 


The  purpose  of  this  study  is  to  predict  the  effect  of 
two  symmetrically  placed  flat  plates  on  the  far  field  of  a 
parabolic  dish  antenna.  The  motivation  for  this  work 
relates  to  the  use  of  dish  antennas  in  environments  where 
unwanted  signals  are  likely  to  be  present.  If  a  simple 
approach  can  be  found  to  predict  the  effect  of  moving  the 
plates  out  from  the  surface  of  the  dish,  then  one  can  on 
command,  move  the  plates  to  an  appropriate  location  so  as 
to  cancel  out  the  effect  of  an  unwanted  signal  by  placing 
nulls  in  the  direction  of  the  unwanted  signal. 

This  study  was  sponsored  by  Daniel  Jacavanco  at  the 
Electromagnetic  Sciences  Division,  Rome  Air  Development 
Center  (RADC) ,  at  Hanscom  Air  Force  Base  Massachusetts  and 
was  part  of  a  continuing  analysis  (Havens,  1983,  Rudisill, 


Background 


Adaptive  antennas  find  their  greatest'  application  in 
situations  where  unwanted  signals  at  the  same  frequency  as 
the  desired  signals  are  being  received.  Such  signals  may 
be  deliberate  or  accidental  but  usually  their  direction  of 
arrival  is  not  along  the  boresight  of  the  receiving 
antenna.  This  means  that  an  antenna  that  can  change  the 
locations  of  its  nulls  in  response  to  the  unwanted  signals 
can  greatly  reduce  the  level  of  the  interfering  signals  and 
yet  receive  the  desired  signal  satisfactorily. 


The  technique  usually  employed  at  present  involves  the 
use  of  a  phased  array  antenna  where  the  phases  of  the 
individual  antenna  elements  can  be  adjusted  to  synthesize 
practically  any  configuration  of  receive  antenna  pattern 
(Johnson  and  Jasik,  1984:9-13).  Selective  nulling  in  the 
receiver  field  pattern  can  be  achieved  through  the  use  of  a 
parabolic  reflector  antenna  having  symmetrically  placed 
phase  modifying  elements.  Standard  reflector  antennas  with 
various  feed  configurations  have  been  thoroughly  analyzed 
both  in  the  near  and  far  field  (Balanis,  1982:604-642; 
Milligan,  1985:220-265;  )  but  symmetrically  placed  flat 
plates  cannot  be  handled  easily  by  the  current  models. 


This  particular  research  topic  has  a  number  of  avenues 


all  of  which  would  prove  to  be  of  unusual  interest. 

However  in  view  of  the  time  constraints  imposed  by  the 
thesis  program  the  scope  will  be  limited  to: 

1.  Scaling  the  reflector  elements  so  that  all 
simplifying  assumptions  are  valid. 

2.  Developing  a  mathematical  model  that 
satisfactorily  describes  the  near  field. 

3.  Correlation  of  computed  to  measured  near  field 
results  and  analysis  of  discrepancies. 

4.  Use  of  the  Fast  Fourier  Transform  to  compute 
far  field  radiation  pattern. 

Approach 

The  approach  in  this  problem  has  been  to  carefully 
analyze  the  desired  physical  setup  to  ensure  that  all 
simplifying  assumptions  are  valid  or  else  to  take  into 
account  their  effect  on  the  accuracy  of  the  final  result. 
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Also  the  area  of  interest  has  been  limited  to  a  region 
extending  thirty  degrees  either  side  of  the  antenna 
boresight.  This  is  where  interfering  signals  have  their 
most  effect  since  the  antenna  being  jammed  still  has 
appreciable  gain  within  this  region. 

Materials  and  Equipment 

RADC  as  the  sponsor  of  this  work  built  a  near  field 
range  so  that  a  model  of  the  total  antenna  system  could  be 
tested  to  obtain  near  field  measurements.  The  results 
obtained  in  the  near  field  range  have  acted  as  a  bench  mark 
to  assess  the  validity  of  the  model  used  for  the  analytical 
derivation  of  the  near  field  pattern. 

Other  Support 

Other  equipment  such  as  computer  support  was  available 
on  site  at  RADC  and  within  the  Air  Force  Institute  of 
Technology  (AFIT) .  These  facilities  proved  to  be  entirely 
satisfactory  except  for  the  fact  that  there  was  no  easy  way 
of  converting  data  obtained  from  the  near  field  range  using 


a  Hewlett  Packard  computer  to  the  MS-DOS  format. 


To  develop  the  far  field  pattern  of  the  parabolic  dish 
antenna  through  the  use  of  a  near  field  range  turns  out  to 
be  no  small  undertaking.  The  dish  to  be  analysed  was  four 
feet  in  diameter  and  the  maximum  height  of  the  measurement 
plane  was  only  52  inches  while  its  width  was  79  inches.  As 
there  was  no  way  to  modify  the  setup  this  limitation  had  to 
be  accepted. 

Unwanted  Reflections 

One  consideration  of  major  importance  in  any  near  field 
measurement  is  the  problem  of  unwanted  reflections  (Newell, 
1985:38).  Use  of  RADAR  absorbent  material  (RAM)  can 
significantly  reduce  the  effects  of  unwanted  reflections 
and  the  chamber  had  been  designed  with  this  criteria  in 
mind.  However  the  vertical  support  that  carried  the 
measurement  probe  was  not  protected  with  RAM.  Since  this 
was  a  piece  of  heavy  guage  steel  approximately  six  inches 
wide,  its  potential  for  reducing  the  accuracy  of  the 
measurements  was  substantial.  A  solution  to  overcome  this 
problem  was  to  offset  the  probe  and  suspend  RAM  material 


from  the  top  of  the  vertical  column  thereby  allowing  the 
measurement  probe  full  unobstructed  vertical  mobility  . 
This  solution  was  quite  successful  but  it  reduced  the 
useable  width  of  the  chamber  by  approximately  eight  inches 


Antenna  to  Measurement  Plane  Separation 


A  further  question  that  was  not  immediately  answered  was 
at  what  distance  should  the  measurement  plane  be  from  the 
antenna.  Initial  measurements  were  carried  out  with  the 
closest  part  of  the  antenna  approximately  two  feet  from  the 
antenna.  As  the  system  was  being  tested  at  3.2  Ghz  this 
was  well  inside  the  near  field  region  with  a  radius  of 
approximately  107  inches  (Balanis,  1982:22).  Analysis 
showed  that  the  pattern  obtained  was  not  providing 
sufficient  information  about  the  antenna  pattern  off  to  the 
sides  (Newell,  1985:39). 


Subsequently  the  antenna  was  moved  so  that  the  nearest 
portion  of  the  antenna  was  within  two  inches  of  the 
measurement  plane.  At  this  distance  the  effect  of  the  feed 
support  became  quite  dramatic  causing  a  drop  of  over  15db 
in  the  near  field  amplitude  and  approximately  50  degrees 
increase  in  the  phase  (See  Figure  2).  From  previous  work 
(Newell,  1985:39),  the  angular  range  over  which  the  far 


field  is  valid  can  be  computed  using  Equation  1. 

Ac  =  tan- 1  (Lx  -  D)  /  (2d)  (1) 

where  Lx  =  Scan  Length  =  79.37  inches 

D  =  Antenna  diameter  =  48  inches 

d  =  Aperture  plane  to  measurement  plane  =  21  inches 
Ac  *  Azimuthal  angle  over  which  far  field  results  are 
valid  =  36.76° 


Near  Field  Errors 


Unfortunately  all  the  measurements  taken  in  the  near 
field  could  not  be  observed  at  the  same  time  and  a  phase 
difference  of  approximately  100  degrees  from  the  top  to  the 
bottom  of  the  projected  dish  was  not  observed  (See  Figure 
3) .  From  Equation  B2  and  the  paragraph  following,  the 
phase  of  the  reflected  field  as  observed  at  the  measurement 
plane  of  a  parabolic  dish  should  have  a  constant  value. 

This  observed  anomoly  could  only  have  been  caused  by  the 
top  of  the  antenna  and  the  measurement  plane  being  closer 
than  the  bottom  due  to  either  misalignment  of  the  dish  or 
tilting  of  the  vertical  support  structure  due  to  the  added 
weight  of  the  RAM. 


During  dish  alignment  a  phase  difference  (across  AA ' 
shown  in  Figure  1)  of  approximately  300  degrees  was  noted. 
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Anomolous  Effects 


Presenting  the  near  field  data  as  a  contour  map  (See 
Figure  4)  showed  an  unusual  ripple.  This  ripple  showed  up 
when  moving  in  a  vertical  direction  through  the  plot.  The 
effect  of  this  high  frequency  ripple  was  to  introduce  an 
apparent  aliasing  effect  (Newell,  1985:68)  when  the  FFT  was 
computed.  The  cause  of  this  effect  could  only  have  been 
due  to  some  anomoly  in  the  sampling  technique  (Reddy, 
1985:9).  A  raster  sampling  technique  was  used  and  the 
program  that  controlled  the  probe  (Appendix  C)  should  have 
moved  the  probe  up  to  the  next  row  in  the  measurement  plane 
and  then  before  moving  on  should  have  taken  the  first 
measurement.  However  apparently  the  probe  moved  one 
position  horizontally  before  taking  a  measurement.  An 
alternative  explanation  is  that  there  was  some  backlash  in 
the  main  worm  drive  that  moved  the  vertical  column  (Newell, 
1985:41).  To  overcome  this  effect,  a  short  routine 
(Appendix  D)  was  inserted  in  the  processing  software  to 
move  every  second  row  of  data  one  position  horizontally 
(See  Figure  5  for  an  example  of  corrected  results).  Once 
this  was  completed  and  the  data  was  processed  the  aliasing 
effect  no  longer  occurred. 
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Large  Dynamic  Range  of  Results 

After  completing  the  measurements  the  processing  of  the 
data  remained  and  during  this  phase  of  the  project  some 
more  interesting  phenomena  were  observed.  One  perplexing 
aspect  related  to  the  dramatic  range  of  the  FFT  results 
over  80db.  After  some  research  a  possible  explanation 
(Newell,  1985:51)  for  this  was  uncovered.  Since  the  near 
field  measurements  were  taken  so  close  to  the  antenna 
evanescent  waves  were  measured.  The  far  field  computed 
using  a  FFT  consists  of  the  sum  of  two  parts  consisting  of 
an  exponential  term  due  to  the  evanescent  waves  and  a 
second  due  to  all  sources  of  measurement  error  (Newell, 
1985:68-69) .  In  Figure  6  an  exponential  amplitude 
distribution  has  been  superimposed  onto  the  far  field 
distribution  to  illustrate  this  point. 

Interpretation  of  Results 

Another  feature  of  interest  was  noted  when  the  data 
points  obtained  from  the  FFT  were  to  be  interpreted.  The 
output  of  the  FFT  was  in  terms  of  Kx/K  and  Ky/K  as  shown  in 
Figure  6  and  not  azmiuth  and  elevation. 


The  values  of  Kx/K  and  Ky/K  range  from  -2.976  to  2. SIS  and 
-2.271  to  2.271  respectively  (See  Appendix  A  page  3).  The 
values  of  Kx/K  and  Ky/K  having  a  magnitude  greater  than  one 
correspond  to  the  transformed  evanescent  waves.  Only  in 
the  regions  where  the  magnitudes  of  Kx/K  and  Ky/K  are  less 
than  one  can  the  output  of  the  FFT  be  used  to  describe  the 
far  field  since  to  relate  the  transformed  data  to  azimuth 
and  elevation,  the  sine  of  the  Kx/K  and  the  Ky/K  values 
respectively,  must  be  taken  (Newell,  1985:6).  Within  this 
range  only  a  somewhat  smaller  angle  can  really  be 
considered  to  be  effective  data  as  found  by  solving 
Equation  1.  To  improve  the  accuracy  of  the  results  the 
size  of  the  measurement  plane  should  be  increased  as  this 
would  effectively  reduce  the  size  of  Kx/K  or  Ky/K  (refer  to 
Equations  A8  and  A9)  and  thereby  lead  to  more  useful  data. 
Zero  filling  can  achieve  similar  results  and  involves 
explicitly  inserting  zeroes  (these  zeros  are  already 
implicit  due  to  the  truncated  area  of  measurement)  in  the 
near  field  data  outside  of  the  areas  where  actual 
measurements  were  taken  (Newell,  1985:28).  Inserting 
zeroes  does  not  imply  any  new  data  and  the  net  effect  is  to 
increase  the  values  of  Lx  or  Ly  which  as  indicated  above 
improves  resolution. 


Data  Sample  Spacing 

Further  gains  could  be  achieved  by  changing  the 
distances  between  samples  to  approximately  half  wavelength 
increments  (Newell,  1985:26).  There  is  a  one-to-one 
mapping  of  points  from  the  near  field  to  the  far  field  and 
using  a  spacing  closer  than  a  half  wavelength  does  not  lead 
to  greater  azimuthal  resolution  in  the  far  field  plot 
(Newell,  1985:28).  Also  it  should  be  noted  that  the  number 
of  sample  points  needed  for  use  in  the  FFT  algorithm  (see 
Appendix  A  for  development  and  Appendix  E  for  the  source 
code  used)  is  an  even  integer  power  of  two. 


Near  Field  Phase  and  Amplitude  Distribution 

Figure  7  shows  the  moveable  plate  as  the  line  segment  RS 
with  GH  being  perpendicular  to  the  tangent  (EGF)  at  G.  As 
the  plate  RS  moves  out  from  the  surface  of  the  parabola  on 
the  axis  GH  it  causes  rays  from  the  focus  f  to  be  reflected 
up  to  the  reference  plane  ABCD.  Figures  8  and  9  show  a  3- 
dimensional  representation  of  the  near  field  phase  for  the 
parabolic  reflector  with  the  plates  at  0  and  225  degrees 
respectively.  The  effect  can  be  clearly  seen  in  the  form 
of  twin  bulges  appearing  symmetrically  either  side  of  the 
bulge  caused  by  the  feed  support.  The  line  segment  BD  in 


Figure  7.  Near  Field  Zones  Due  to  Moveable  Plate 

figure  7  denotes  the  region  of  illumination  due  to  the 
moveable  plate  which  moves  to  the  right  as  the  plate  moves 
out  from  the  surface  of  the  parabola.  Line  segment  AB 
represents  the  region  on  the  measurement  plane  that  is 
shadowed  by  the  flat  plate.  As  the  flat  plate  moves  out 


from  the  parabolic  surface  AC  increases  in  length  and 
region  AB  is  not  illuminated  by  either  rays  reflected 
directly  off  the  parabolic  dish  or  rays  reflected  off  the 
moveable  plate  (See  Appendix  E  for  a  program  that  computes 
the  co-ordinates  of  the  various  regions  of  this  second 
model  in  the  near  field  measurement  plane) .  Rays 
reflecting  off  the  plate  illuminate  the  region  BD,  but 
region  CD  is  also  illuminated  by  rays  reflecting  off  the 
dish. 

Figures  8  to  16  show  the  distribution  of  energy  in  the 
the  near  field  while  Figures  17  to  25  show  the  phase 
distribution  as  the  plates  are  moved  out.  These 
distributions  are  those  recorded  as  the  measurement  probe 
travelled  over  the  slice  AA'  shown  in  Figure  1.  The  phase 
distributions  show  the  twin  bulges  as  expected  and  as  the 
plates  move  out  the  phase  disturbance  increases  until  225° 
is  reached.  After  this  it  decreases  and  at  360°  the  bulges 
dissapear  entirely. 

Effect  of  Moving  Plates 

If  diffraction  is  ignored,  then  as  the  plate  (refer  to 
Figure  7)  moves  out  (in  1/16  wavelength  increments)  from 
the  dish  surface,  the  transit  time  required  for  the  rays 
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emaninating  from  the  focus  to  travel  to  the  measurement 
plane  is  reduced.  When  the  plate  is  at  1/2  wavelength  from 
the  dish,  a  ray  emaninating  from  the  focus  striking  the 
plate  and  reaching  the  measurement  plane,  has  travelled  a 
shorter  distance.  This  path  is  shortened  by  approximately 
one  wavelength  over  a  ray  that  struck  the  dish  and 
reflected  to  the  measurement  plane  directly.  Thus  these  two 
rays  are  approximately  back  in  phase.  This  effect  can  be 
seen  quite  clearly  by  comparing  Figures  17  and  25  which 
show  the  near  field  phase  distribution  at  0  and  360  degrees 
respectively  to  be  of  fundamentally  the  same  shape. 


Computed  Near  Field  Phase  Distribution 


The  Near  Field  Phase  Distribution  Program  (Appendix  E) 
calculates  the  near  field  phase  distribution  for  one  plate 
since  the  results  on  the  other  side  of  the  feed  blockage 
should  be  a  mirror  image.  This  is  an  admittedly  simple 
approach  which  relies  entirely  on  ray  tracing  techniques 
and  vectorial  summation  of  rays  reflected  directly  off  the 
dish  and  those  rays  reflected  off  the  moveable  plate  to 
model  a  very  complicated  problem.  Modelling  only  a  small 
portion  of  the  near  field  is  considered  reasonable  since 
one  of  the  stated  aims  of  this  study  was  to  develop  a 
predictive  ability  regarding  the  far  field  behaviour  as  a 


function  of  plate  movement.  Apart  from  being  a  two 
dimensional  model  no  attempt  is  made  to  take  account  of 
diffraction  effects  which  certainly  would  be  important  in 
view  of  the  relatively  small  size  of  the  moveable  plates. 
Also  the  feed  blockage  causes  a  severe  effect  on  the  near 
field  phase  distribution  and  this  has  not  been  included. 

Considering  the  results  obtained  there  is  an  apparent 
correlation  between  the  measured  and  computed  near  field 
phase  distributions.  This  qualitative  correlation  seems 
strongest  in  the  regions  where  the  horn  blockage  and 
diffraction  effects  would  be  expected  to  be  smallest  (ie. 
nearest  to  the  outside  edge  of  the  dish) .  In  most  cases 
the  abrupt  transition  from  the  larger  phase,  attributed  to 
the  moveable  plate,  back  to  the  smaller  phase,  due  to 
reflection  directly  off  the  surface  of  the  dish,  in  the 
computed  results  coincides  reasonably  in  position  and 
magnitude  with  the  transition  apparent  in  the  measured 
plots.  These  results  give  some  confidence  to  the  concept 
of  this  simplified  model.  Adding  to  this  model  those  rays 
which  are  diffracted  by  the  edge  of  the  plate  may  further 
enhance  the  accuracy  of  the  predicted  near  field 
distribution . 

A  finite  identifiable  number  of  rays  diffract  directly 


off  the  moveable  plate  and  could  strike  a  single  point 
along  AA'  (Figure  1) .  There  are  two  sources  for  these 
diffracted  rays.  The  first  source  is  rays  that  strike  the 
plate  edge  with  normal  incidence  and  the  second  is  rays 
that  strike  the  plate  edge  with  oblique  incidence.  These 
rays  with  appropriate  diffraction  coefficients  would  serve 
to  enhance  this  simple  model. 
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Figure  19.  N/Field  Phase  Distribution  -  Plates  at  0° 
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Figure  21.  N/Field  Phase  Distribution  -  Plates  at  90 


PM. 


ly  *-l  I. 


Figure  24.  N/Field  Phase  Distribution  -  Plates  at  225° 
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Figure  26.  N/Field  Phase  Distribution  -  Plates  at  315® 


IV.  Comparison  of  Measured  to  Computed  Results 


The  near  field  results  shown  in  the  previous  chapter 
were  processed  using  the  fast  fourier  transform  program 
shown  in  Appendix  D  and  sample  amplitude  distributions  are 
provided  by  figures  28  to  32  inclusive.  Plots  of  the 
corresponding  measured  data  is  shown  in  figures  33  to  37. 

Comparison  of  Side  Lobe  to  Main  Lobe 

Comparison  of  the  results  is  facilitated  by  considering 
the  amplitude  of  the  first  side  lobe  relative  to  the  main 
lobe.  In  each  case  reasonable  correlation  of  the  measured 
to  computed  results  is  obtained.  (See  Table  1).  Where 
discrepencies  occur  it  is  an  indication  of  two  problem 
areas.  The  first  is  the  size  of  the  measurement  plane 
which  could  be  extended  by  zero  filling  to  provide  enhanced 
resolution  in  the  far  field  (Newell,  1985:28).  Also  the 
taking  of  accurate  measurements  in  the  near  or  far  field  is 
made  very  difficult  by  the  large  number  of  systematic 
sources  of  error.  These  error  sources  are  (Newell, 
1985:67-85)  as  follows: 


1. 


Multiple  reflections. 


2.  Probe  positioning  uncertainties. 

3.  System  amplitude  and  phase  nonlinearities  (receiver 
and  attenuator) . 

4.  Impedence  changes  due  to  probe  and  receiver  cable 
movement . 

5.  Aliasing  if  no  band  limit  exists  on  the  FFT 
results  {band  limit  is  inversely  proportional  to 
the  size  of  the  near  field  sample  step) . 

6.  Aliasing  if  near  field  samples  are  spaced  too 
widely. 

7.  Area  truncation  (ie  sampling  over  a  limited  area  in 
the  near  field) . 

Null  Movement 

Also  of  interest  is  how  much  does  the  first  null  move  as 
the  plates  move  out  from  the  surface  of  the  dish.  The 
computed  results  do  not  provide  sufficient  resolution  to 
allow  any  meaningful  conclusions  to  be  drawn.  However  the 
measured  results  do  indicate  that  there  was  some  slight 
movement  of  the  null.  (See  Table  2). 
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Figure  31.  Computed  Far  Field  Pattern 
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V.  Recommendations  and  Suaaestions  for  Further  Work 


Areas  Worthy  of  Further  Investiaation 


During  the  course  of  the  investigations  into  the  effects 
of  moveable  flat  plates  on  the  surface  of  a  parabolic  dish 
antenna  some  areas  where  further  investigation  would  be 
warranted  were  discovered.  Due  to  time  considerations 
these  potentially  interesting  areas  could  not  be 
investigated.  Aspects  that  may  well  prove  worthy  of 
further  work  include: 


Analyzing  whether  parasitic  fields  are  set  up 
between  the  reflector  disks  and  the  main 
reflector  and  if  so  how  they  influence  the 
final  results. 


2.  Considering  the  effects  that  different  field 
distributions  in  the  feed  horn  have  on  the 
near  field  pattern. 


3.  Introducing  the  effects  that  random  variations 
in  the  surface  of  the  main  reflector  will 
cause  to  the  near  field. 


4. 


Reducing  the  size  of  the  parabolic  dish  so 
that  the  region  of  interest  is  not  on  the 
extreme  of  the  useable  data  from  the  FFT. 


5.  Using  disks  that  can  be  tilted  with  respect  to 
their  support  shaft  so  that  the  angular 
location  of  the  null  might  be  moved. 

6.  Increasing  the  size  of  the  disks  to  determine 
if  further  enhancement  of  the  null  may  be 
possible  or  if  degradation  of  the  main  lobe 
would  be  the  major  effect. 

Implementation  of  Moment  Method  Technique 

In  all  moment  method  codes  there  are  essentially  three 
main  problems  that  must  be  addressed  before  valid  results 
can  be  obtained.  Firstly  the  structure  must  be  gridded  up 
so  that  no  segment  has  a  side  longer  than  1/4  wavelength 
and  the  end  point  co-ordinates  for  all  of  these  separate 
segments  must  be  fed  to  the  moment  method  program.  For 
anything  more  than  very  simple  structures  this  procedure 
must  be  handled  using  a  computer  program.  Secondly  since 
only  straight  lines  allowed  by  the  code  the  most 
appropriate  method  of  approximating  curved  structures  must 


be  considered.  Finally  a  model  of  the  antenna  feed 
structure  that  is  able  to  approximate  the  salient  features 
of  the  feed  distibution  is  required.  This  last  requirement 
can  cause  some  difficulty  since  the  model  should  account 
for  the  feed  polarization  as  well  as  the  overall  pattern. 

Whether  a  simplified  model  can  be  derived  that 
approximately  accounts  for  all  observed  effects  or  whether 
the  method  of  moments  is  the  only  technique  that  can 
provide  the  required  accuracy  is  unclear  at  the  moment. 


Appendix  A:  Development  of  the  FFT  Algorithm 
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where  M  =  2  and  N  =  2  where  M  and  N  are  the  number 
of  rows  and  columns  in  the  reference  plane  respectively. 
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then 


where  p  =  0,1,2, . ,M-1  and  q  =  0,1,2, 
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To  put  this  into  the  standard  form  of  an  FFT  define  W 


and  W  as  follows: 
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however  in  the  discrete  Fourier  integral  above  (Eq  A7) 


we  find  that: 
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substituting  in  the  values  for  A  and  A  produces 
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substituting  these  values  for  W  an d  W  into  the  discrete 

M  Pi 


Fourier  integral  given  above  (Eg  A7)  results  in  the  following: 
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The  implementation  of  Equation  A10  is  shown  in  Appendix  E. 


Equations  A8  and  A9  permit  the  determination  of  the  maximum 


values  of  A  and  A 
X  y 


p  =  M— 1  =  127 

ma  x 


L  =  1.008  metres 

X 


q  =  N— 1  =  63 
ma  x 


L  =  . 66  metres 


k  =  2«/A 


A.  =  c/f 


c  =  3x10  metres/sec 


f  =  3.2x10  Hz 


Substituting  the  above  values  into  Equations  A8  and  A9  the 


following  values  for  k  /k  and  k  /k  are  obtained 

x  y 


k  /k  e  [-2.976,  2.976]  and  k  /k  e  [-2.271,  2.271] 

x  y 


Appendix  B:  Development  of  an  Analytical  Model 


Consider  the  situation  firstly  from  the  point  of  view  of 
tracing  rays  from  the  focus  of  the  parabolic  dish  antenna  to  the 
surface  of  the  dish. 


Suppose  that  the  focus  is  at  the  origin  and  the  distance  to 
the  directrix  of  the  parabola  is  d.  Since  this  model  is  to  be 
compared  to  the  experimentally  derived  results  let  us  consider 
that  there  is  a  measurement  plane  a  distance  ft-  away  from  the 
origin  and  parallel  to  the  x-axis  as  shown  in  Figure  38. 

We  can  now  write  the  equation  of  the  parabola  in  polar  form  as 
follows : 


r  -  g{ 0) 


d 

1  -  sin  Q 


(Bl) 


where  r  is  the  distance  from  the  focus  to  the  surface  of  the 
parabolic  surface  and  0  is  the  angle  measured  clockwise  from 
the  x-axis. 

Now  to  define  the  distance  of  the  reference  plane  from  the 
surface  of  the  parabola  we  need  to  first  note  that  the  distance 
from  the  x-axis  to  the  parabolic  surface  is  -r  sin  0  therefore 


the  distance  to  the  reference  plane  is 


ft  +  ( -r  sin 


therefore  the  distance  travelled  by  any  ray  from  the  focus  to  the 
measurement  plane  is 


D  =  r  +  h  -  r  sin  B 


D=h+r  ( 1  -  sin  0) 


Using  Equation  (1)  here,  we  obtain. 


D  =  h 


d  (1 

1  -  sin  B 


sin  B ) 


D  =  h  +  d 
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D  is  a  constant  value  independent  of  B  which  was  expected 
since  a  wave  of  constant  phase  should  be  reflected  from  the 
surface  of  a  parabolic  surface. 


The  next  step  is  to  find  the  gradient  of  a  tangent  to  the 
surface  of  the  parabola  at  any  point.  Note  that  this  will  be  the 
gradient  of  the  moveable  flat  plate  which  moves  up  from  the 
surface  of  the  parabola  since  it  is  mounted  on  a  shaft  which  is 
perpendicular  to  the  surface  at  this  point.  (See  Figure  38) 


Suppose  that  the  angle  from  positive  direction  of  the  x-axis 


to  the  general  point  (x,y)  is  B  then 
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Considering  the  situation  depicted  in  Figure  (40). 

The  next  step  is  to  find  the  equation  defining  the 
line  of  which  the  flat  plate  is  a  line  segment.  Let  (  <  ,  0  ) 
be  the  polar  co-ordinates  of  a  general  point  on  the  line. 

The  line  segment  can  be  defined  to  exist  between  9  and  &2  ie. 


0  e  iei  ,  e2] 

From  equation  (B3)  we  have  the  gradient  of  the  line  segment 


y  -  K  sin  0 
x  -  K  cos  0 


=  tan  a 


( B8 ) 


into  which  the  equations  in  polar  co-ordinates  for  x  and  y 

can  be  substituted  as  follows: 

r  sin  0  -  K  sin  <t>  _  sin  a 

r  cos  0  -  <  cos  0  -  cos  a 

cos  a  (r~  sin  9  -  K  sin  0)  =  sin  a  (r  cos  9  -  t,  cos  0) 


r  (cosci  sin0  -  sina  cos^) 


K  (cosa  sin0  -  sina  cos0) 


r  sin  (a  -  0)  =  t,  sin  (a  - 
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Also,  from  (B8)  r 
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Considering  next  a  constant  phase  ray  as  it  scatters  from  the 


parabola  described  by  equation  (Bl) 
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where  r  =  g($)  (cos  9,  sin  &> 


Next,  consider  a  constant  phase  front  ray  as  it  scatters 
from  the  line  segment  described  by  equation  (B9). 


h(0)  = 


C  sin  ( o  -  0) 
sin  ( ct  -  0 ) 


r  ( t )  = 
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where  s  -  h ( @ )  (cos  sin  £> 

and  3  is  defined  in  equation  (B6)  as  3  =  20  -  <t> 


However  if  the  line  segment  is  defined  to  lie  between  9  s 
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then  we  have 
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where  the  distance  from  the  x-axis  to  the  reference  plane  is 


d  =  Im  {  h  ( 3 )  e  J  ^  +  (ct  -  h(0))eJJ>} 


Itf. 


=  h(£)  sin  0  +  (ct  -  h(£))  sin  (2 ct  -  0) 


and  the  path  length  from  the  focus  to  the  reference  plane  is 


d  -  h(0)  sin  0  .  , ^ 

ct  =  - : - r— - c-t—  +  h(©J 

sm  (2Ci  -  0) 


The  next  step  is  to  find  the  boundaries  of  the  regions  indicated 
in  Figure  43.  Noting  that  the  center  of  the  plate  is  defined  by 
T(£,0),  its  cartesian  co-ordinates  are 
x  co-ordinate  =  K  cos  0 
y  co-ordinate  =  K  sin  0 

now  the  co-ordinates  of  R  and  S  can  be  written  in  cartesian  co¬ 


ordinates  as  follows: 


R(<  cos  0  -  k/2  cos  ct,  K  sin  0  -  k/2  sin  a.) 

S(K  cos  0  +  k/2  cos  a,  K  sin  0  +  k/2  sin  a) 

where  k  is  the  diameter  of  the  moveable  plate  and  ct  is  the  gradient 


of  the  plate. 


Now  finding  the  co-ordinates  of  B  and  D: 


K  sin  0  -  k/2  sin  a 


let  Yj  =  (  sin  0  +  k/2  sin  ct  and  y2  = 

where  y  and  y2  are  the  y  co-ordinates  of  S  and  R  respectively 

D  ( ^  cos  0  +  k/2  cos  ct  +  (yi  +  D)/tan  3  ,  D) 

B(C  cos  0  -  k/2  cos  cl  +  ( y  ^  +  D)/tan  #2,  D) 

A(£  cos  0  -  k/2  cos  ct,  D) 

C(r  cos  0  ,  D  )  where  r  =  d/Cl  -  sin  0)  and  3  and 

3  are  the  angles  that  rays  bouncing  from  the  extreme  edges  at  S 

R  respectively  make  with  the  positive  x  direction. 

The  point  T  is  defined  by  its  distance  above  the  surface  of  the 
parabola  (Bl)  therefore  T  needs  to  be  available  in  terms  of  this 
variable . 

The  gradient  of  the  tangent  is  cos  0Q/(1  -  sin  0  )  therefore  the 

gradient  of  the  shaft  to  which  the  plate  is  attached  will  be 

7  =  tan  1(  -(1  -  sin  0Q)/cos  <?Q) 
where  7  is  the  angle  that  the  perpendicular  makes  with  the  x-axis 

therefore  the  co-ordinates  of  T  are  as  follows: 

T(r  cos  0-1  sin  (/  -  90),  -r  sin  9  -  l) 


Appendix  C 


:  Near  Field  Measurement  Control  Program 


'HIS  =ROQRAM  WAS  USED  '0  CONTROL  ’’HE  4EAA  FIELD  EXPERIMENT  CONDUCTE! 
-OF.  GLEN  'HC.RFE  IN  THE  SL'TMEF  ZF  :=5'. 

THE  r  F.CGFj-i’!  : 

mj  C3NTr\CL^  'HE  'HO  DIM  SCAN  .  a  ,  "  )  OF  i-N  ArERA'_RE  -FORE. 

B>  READ-  -ND  =TOFE=  DATA  FROM  'HE  SCIENTIFIC  AT— .-iNTA  *78  j 
FHHSE  AMPLITUDE  RECEIVER  .  CONTROL  OMIT. 

C,  CONTROLS  THE  MOVEMENT  CF  TEC  DISC  LOCATED  CM  A  FOUR  FCCT 
DISH. 


OPTION  BhSE  I 

DIM  A>t,p  UC3)  .Phase':  109) 

PRINTER  IS  I 

I NT  EGER;  r  t.5  ,  No  ,  Vdos  ,  Ho  os  ,  Hstpc  .  Vstao  .  Di  r  -  .  Di  r  v  .  Hstp  .  »' stp  .  Mot  :r  ,  K 

'The  next  rsw  lines  are  uses  to  se~  of'  the  run 


PRINT  "THE  FFQE-E  MOVES  HORIZONTALLY  -ROM  ="CS1'LEFT1  TO  PCS  103  <  RIGHT ) 
INPUT  "INPUT  THE  PRESENT  HOP  LOCATION  "..-Iota 

PRINT  "'HE  RC05E  MOVES  /E.RT  I  0Ai_L  <  FROM  ROS  lv LOWEST.  TO  PCS  o4  i  H I GHEST ) 
INPUT  "INPUT  'HE  PRESENT  VERT  LOCATION  ",Mor& 

INPUT  "INPUT  THE  PRESENT  DISC  POSIT  ION" . Dsc 
INPUT  "DO  YOU  WANT  TO  MOVE  DISC  BEFORE  THE  RUN  <Y/N)" ,YnS 


IF  Ynf  =  "  M"’  THEN  210 
CALL  Dsc.ncv  (  Dsc  > 

INPUT  "DO  YCU  WISH  TO  SET  POWER  AT  A  SPECIFIC  LOCATION 
IF  AS- "N"  THEN  550 
IF  AS  .  "  THEN  210 


ON  KEY  5  LABEL 
ON  KEY  6  LABEL 
ON  KEY  0  LAE EL 
ON  Y EY  7  LABEL 
CM  KEY  1  LABEL 
GOTO  240 


" F ROBE  LEFT"  GOTO  TOO 
"PROBE  DOWN"  GOTO  420 
"RET.  TO  MAIN"  GOTO  550 
"PROBE  RIGHT"- GOTO  040 
"PROBE  UF"  GOTO  030 


(Y/N) " .AS 


* 


A  =  4 


Sto  =  4R>I 
Mota=Mota- V 
GOTO  420: 

A=  0 

3tp=4R : 
Mota=Mci: a-  i 
GOTO  420 


S  t  s  =  7  4  2 

M 0-7)=  More  i-  1 


V 


FOR  ,L=1  TO  Stp 

OUTF'-JT  sGcic  USING  "4.  W"  ;  Motor 
CL,_F'UT  55c  1  =:■:-■ 

WAIT  . 0023 
NEXT  K 

PRINT  '  HCR  *  <-CNThL  _CCAT  ION  —  * ; Mo t a ■ "VERT  I  CAL  *_.jCATI3 
LOCAL  TOT 
GC^O  240 

INPUT  "DC  YOU  WISH  TO  PREPOSITION  THE  PROBE v y / N Af 
IF  Af="N"  THEN  590 
IF  A*-.  .  "Y"  THEN  550 


ION  = " : Mots 


580 

CALL  .1ct.TicviTtota.Motj 

590 

vpos=e4 

600 

D: r v=L 

610 

Vstpp= 1 

620 

Di rn=T 

630 

Hoos=128 

1  640 

Hstpp— 1 

I  650 

E=Mota 

1  660 

IF  Motb=o4  THEN 

1  670 

Di rv= 1 

680 

Vpos=l 

690 

Vstpp=-1 

700 

END  IF 

710 

IF  Motb=l  THEN 

720 

Di  rv=2 

730 

Vpos=64 

740 

Vstpp=l 

750 

END  IF 

760 

F=Motb 

770 

ON  KEY  2  LABEL  "END  FROG" 

GOTO  1630 

.  780 

FQR  Posb-F  TO  Vdos  STEP  Vstpp 

790 

IF  Mota=128  THEN  Dirh=4 

800 

IF  Mota=12S  THEN  Hpos=l 

810 

IF  Mota=128  THEN  Hstpp= 

-1 

820 

IF  Mota= 1  THEN  Dirh=T 

830 

IF  Mota=l  THEN  Hpos=128 

840 

IF  Mota=l  THEN  Hstpp=l 

850 

FOR  Posa=E  TO  Hpos  STEP 

Hstpp 

860 

Mota=Posa 

870 

Motb=Poso 

880 

PRINT  "HORIZONTAL  =": 

Mot a, "VERT 

881 

1 

882 

! THE  NEXT  SET  OF  CODE 

READS  THE 

883 

1 

890 

ASSIGN  '©Gpio  TO  12: FORMAT  OFF  ! 

900 

ASSIGN  <SR  TO  70 3 

910 

REMOTE  SR 

920 

OUTPUT  @R:"W1" 

930 

OUTPUT  SR: "G8To5" 

940 

CLEAR 

950 

WAIT  .1 

960 

i  p  :  T  F  I G  G  E  F  ‘3  R 

9  7  <  j 

St  at:  I  =  tat=9P0LL  i  ?'• 

980 

Ir  BI_.Istat,";=:  AND 

B  .  _ '  1st jt 

990 

I-  BIT  .  !  st  Jt .  '  =:  'HE 

3  L,  '• 

1  00(! 

I  -  B  I  ~  I  at  at .  5  ;  THE; 

l  Bt  at 

1  0  i  •! 

s. ;  E-  iP  :  A  .  B  .  2  .  wc  .  B  p 

i  0  Z‘. 

H'ImC  ~  )  —A 

i  04 1' 

F  riase i F  osa  > =Ac 

r  OUND  '* 


i  OoO  Sod :  .3 

1070  5er  =SF CLL  i  705; 

1030  PR  I  NT  1  'ESS  I  7ER  15  WCftT  ING  I.'(CCftftc2T 

ioro  local  g 

1100  PRINT  ’ F F ESS  CONTINUE  ONCE  EPRCR  15 

lilO  BEEF 

1120  PAUSE 


FFCP  NUMB 


1130 
1140 
1  ISO 
1151 
1132 
1 153 
1 160 
1170 
1180 
1190 
1200 
1210 
1220 

1230 

1231 

1232 

1233 
1240 
1230 
1260 
1270 
1280 
1290 
1300 
1310 
1320 


CLEAR  2ft 
GOTO  Tq 
Okay:  ! 


'NEXT  MOVE  THE  PROBE  HCEI ZONTALLY 

Motor  =  2'"(Dirh-l) 

FOR  k=l  TO  491 

OUTPUT  @Gpio  USING  " 4, W"; Motor 
OUTPUT  SGpio;0 
WAIT  .0025 
NEXT  K 
WAIT  .3 
NEXT  Posa 
» 

‘NEXT  STOP  THE  DATA  IN  THE  PROPER  MASS  STORAGE  DEVICE 

i 

IF  D=c=360  THEN 

CALL  Datat 1 (Amp ( *) , Phase \ * ) , 129 . Motb , Dsc  > 

END  IF 

IF  Dsc=225  THEN 

CALL  Datat 2 (Amp (■*  >  , Phase (* )  ,  123 , Motb , Dsc ) 

END  IF 

IF  Dsc=3 1 5  THEN 

CALL  DataP3 (Amp (*) , Phase ( *> , 1 28 . Motb , Dsc ) 

END  IF 


1330  Motor  =  2's  (Dirv-1 ) 

1331 

1332  'NOW  MOVE  THE  PROBE  VERTICALLY 

1340  FOR  K= 1  TO  742 

1330  OUTPUT  <sGp  i  o  USING  "».W";Motor 

1360  OUTPUT  illGpiosO 

1370  WAIT  .0025 

1380  NEXT  K 

1390  E=Mota 

1400  NEXT  F'osb 

1410  IF  Dsc=360 'THEN  1610 

1420  Dsc=Csc+45 

1430  Dsca=l 

1440  Dscb=4 

1450  ASSIGN  sGoio  TQ  15: FORMAT  OFF 
1460  Mot or -2  < Dsca- 1 ) 

1470  FOR  k=t  TO  3c3 

1 430  OUTPUT  iGoio  USING  ’’ tt  ,  W"  :  Motor 

1490  OUTPUT  OGp  i  o :  0 

1300  WAIT  .  0023 

15 lv  NEXT  ■ 

1320  Motor=Z  'Dscb-L' 

.  530  FOR  (  —  1  TO  3ijo 

1 340  OUTPUT  aGo i o  USING  "4,W": Motor 

. 350  OUTPUT  aGo:o:0 
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l  a  • »  «  t  i  '|  . 


Ll'(U  U  *  »  !■»  |  •  1^| 


WA  l  T  ,002b 
NEa T  >• 

ASSIGN  aGj  i  :■  ’’"G  :  2; FORMAT  DF- 
FRINT  '‘^h^  OLRRENT  DISC  FG3  I T I ON  IS"  •  D  -  C 
GC7C  joO  1  ONE  MCF.E  TIME 
r  ft  IN7  "  OhTr.  :s  CC«F'.E'E  * 

BEEF  DO l  M-! ,  .  G 
END 

SUB  Datir 1  Amo  >  *  1  . Phase ■ * ;  , Not s .  Motb .  Dsc  J 

'THIS  SUBPROGRAM  5E7S  UP  A  DATA  FILE  CONS IS7 ING  OF  THE  AMPLITUDE  AND  * 
!  PHASE  DATA  RECORDED  IN  THE  MAIN  PROGRAM  * 


Dsc)?."  :  INTERNAL,  4,0"  ,Npts,16 


CREATE  BOAT  "  D_  '  ‘-.VAL*  (Motb  )  ?<"  _"!<VALi  ( Dsc )  S."  :  INTERNAL  ,  4 , 0"  ,Npts,  1 
ASSIGN  SP  TO  "D  _"?<VALS  t  Motb  >  _“?<VAL.p  <Dsc )  -V' :  INTERNAL  ,4,0" 

OUTPUT  aP : Amp (*)  , Phase (* ) 

ASSIGN  aP  TO  * 

PRINT  "DATA_" ; Motb ; "HAS  BEEN  CREATED" 

SUBEND 

SUB  Motmov  (Mota , Motb ) 

ASSIGN  'SGoio  TO  12;  FORMAT  OFF 

INPUT  "INPUT  THE  NEW  HORIZONTAL  LOCATION  FOR  THE  PROBE", Hnew 
Hstp= 1 
Di rh  =  D 

INPUT  "INPUT  THE  REQUIRED  VERTICAL  LOCATION  FOR  THE  PROBE", Vnew 

Vstp=l 

Di rv=2 

IF  Mota-Hnew.-O  THEN  Hstp=-i 
IF  Mota-hnew ,0  THEN  Dirn=4 
IF  MotL-Vnew 0  THEN  Vstp=-1 
IF  Motb-Vnew-0  THEN  Dirv=l 
IF  Mota-Hnew=0  THEN  1970 
FOR  Mota-Mota  TO  Hnew  STEP  Hstp 

PRINT  "HORIZONTAL ;Mota, " VERTICAL= " : Motb 
Motor  =  2'''  (Dirh-1 ) 

FOR  K= 1  TO  491  « 

OUTPUT  aGoio  USING  "#,W":Motor  , 

OUTPUT  aGpio;0 
WA IT  , 0025 
NEXT  K 
NEXT  Mot a 

IF  Motb-Vnew=Oi  THEN  2070 

FOR  Motb=Motb  TO  Vnew  STEP  Vstp 

PRINT  "HCR I ZONTAL= " ; Mot a , "VERTICAL " ; Motb 
Motor =2 "  (Dirv- 1 / 

FOR  K=f  TO  742 

CUTFUT  agolo  USING  "  4,  W“;  Mot  or  * 

OUTPUT  aGpiosO 
WAIT  .0025 
NEXT  t 
NE<T  Morb 
SUBEND 

'THIS  SUBFCCGRh-M  IS  '.'SED  TO  MOVE  THE  TWG  DISC  LOCATED  ON  THE  FOUR 
'  FCOr  DISH 

SLB  Dst tis  /  '  D  =  c  > 

-(til  30  ?:Go  l  j  TO  l  5 ;  F  DF.MA  OF” 

CM  (!•  I  L/-B  EL  "LI  EC  IUT"  3  J  TC  ImO 
ON  r-Er  ’  i_ABE:_  "DISC  IN"  GOTO  213'' 
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p-.T-s -wn  w-y_w  w^rr  TX7~T^r^  r 


2 1 20 

3N  ■  S':  4  LhBEL  "  A1 

2 1  CO 

■SOT  3 

2 1 40 

Dsc  a— 1 

2150 

Qscp=“ 

2 1 60 

8bc  =  -kAC 

2170 

GOTO  221-1 

2  1 80 

Bbo=-AS 

2190 

Dsca=2 

2200 

Dscp=3 

2210 

Mot  or  =  2  < Dsca- 1 > 

2220 

FOR  T=1  TO  56 a 

2230 

OUTPUT  aGoiJ  US 

2240 

OUTPUT  SGpio;0 

2250 

WAIT  .0025 

2260 

NEXT  K 

2270 

Motor  =  2''(Dscb-l) 

2280 

FOR  K= 1  TO  5o8 

2290 

OUTPUT  iSGoio  US 

2300 

OUTF  UT  i§Gp  l  o ;  0 

2310 

WAIT  .0025 

2320 

NEXT  K 

2o>jO 

Dsc  =  Qsc-*-Bbo 

2340 

PRINT  "THE  CURREN' 

2350 

GOTO  2100 

2360 

ASSIGN  AGpio  TO  1. 

2370 

SUBEND 

2371 

'THE  NEXT  TWO  SUBPROGRAMS  ARE  USED  TO  STORE  DATA  IN  DIFFERENT  MSI 

SUB  Da  tat  0  l  Amo  t  *)  ,  Phase  <•*  >  ,  Npts  .Motb .  Dsc  ) 

CREATE  BOAT  "D_"^VALf  (Motb  >  •*< "  _"?<VAL-f  <  Dsc  )  ?<" :  INTERNAL,  4, 1 " ,Npts, 16 
ASSIGN  <iP  TO  "D_“  !<VALf  < Motb  >  V_"S<VAL#  ( Dsc )  S<”  :  INTERNAL ,4,1" 

OUTPUT  SF's  Amp  (*>  .Phase  <*> 

ASSIGN  '«F'  TO  * 

PRINT  "DATA_" ; Motb; "HAS  BEEN  CREATED" 

SUBEND  4 

SUB  Da tat  3  <  Amp  <*>  , Phase  <*) , Npts ,Motb ,  Dsc)  , 

CREATE  BDAT  "D_".!<VAL$  (Motb  )  ?<"_"S<VAL£  (Dsc  )  V  :  C330, 700"  ,  Npts,  16 
ASSIGN  'SP  TO  " D_"$<VAL-J  (Motb >  J<"_"!<VAL-P  (Dsc >•!■<" :  C3S0 , 700“ 

OUTFUT  iSP s  Amp  ( * )  ,  Phase  v  *i 
ASSIGN  <srp  TO  * 

PRINT  "DATA_" ; Motb; "HAS  BEEN  CREATED" 

SUBEND 


*’«*■  V  -■  *y  ’Jr's 

■a**.  ** 


.Lvv,y 


Appendix  D.  Near  Field  Error  Correction  Program 


PRINT  "***********  NFPHAD1 .  BAS  ********************  » 
PRINT  "THIS  PROGRAM  TAKES  THE  SEPARATE  NEAR  FIELD  DATA  " 
PRINT  "FILES  READS  IN  THE  DATA  AND  OUTPUTS  THE  DATA  TO  A  " 
PRINT  "RANDOM  ACCESS  FILE.  THIS  RANDOM  ACCESS  FILE  CAN 
THEN"  PRINT  "BE  USED  THE  PLOTTING  ROUTINES  TO  CREATE 
TOPOGRAPHICAL"  PRINT  "MODELS  AND  SURFACE  REPRESENTATIONS  OF 
THE  NEAR  FIELD"  PRINT  "DATA.  HOW  MANY  DEGREES  ARE  THE 
PLATES  ABOVE  THE  DISH";  INPUT  P% 

PRINT  "YOU  MUST  ASSIGN  A  VALUE  OF  EITHER  O  OR  1  TO  K  " ; 
PRINT  "DEPENDING  ON  WHETHER  YOU  WANT  EVERY  SECOND  ROW  OF  " ; 
PRINT  "NEAR  FIELD  DATA  TO  BE  SHIFTED  TO  THE  LEFT  OR  RIGHT"; 
PRINT  "WHAT  VALUE  WOULD  YOU  LIKE  TO  ASSIGN  TO  K  (0/1)  "; 
INPUT  K% 

INPUT  "HOW  MANY  FILES  OF  INPUT  DATA  ",M% 

MM%=M%-1 

mu%  =  cint (log (m%) /log  (2) ) 

INPUT  "HOW  MANY  DATA  ELEMENTS  IN  EACH  INPUT  FILE  "  ,  N% 
NN%=N%-1 

nu%  =  cint (log (n%) /log (2) ) 

DIM  AMOP! (0:MM%,  0 : NN% ) , PHOP ! ( 0 : MM% ,  0:NN%) 

DIM  AMP! (0:NN%) ,  PH!(0:NN%),  PHASE !( 0 : NN% ) 

PRINT  "***************  WORKING  *****************" 

ALO !  =  9.9E+30 
AHI !  =  -ZLO 
BLO!  =  9.9E+30 
BHI!  =  -ZLO 
IF  P%  <  10  THEN 

P$  =  RIGHTS (STR$ (P%) , 1) 

ELSEIF  P%  >  9  AND  P%  <  100  THEN 
P$  =  RIGHTS (STR$ (P%) , 2) 

ELSE 

P$  =  RIGHTS (STR$ (P%) , 3) 

END  IF 

FOR  R%  =  0  TO  MM% 

IF  R%  <  9  THEN 

B$  =  RIGHTS (STRS (R%+1) , 1) 

ELSE 

B$  =  RIGHTS (STRS (R%+1) , 2) 

END  IF 

c$  =  "A: "+P$+"DAT"+B$+" .BAS" 
open  c$  for  input  as  #1 
K%  =  K%  +  1 
C%  =  0 

DO  UNTIL  EOF ( 1 ) 
input  #1,  AMP! (C%) ,PH! (C%) 

IF  AMP ! ( C% )  <  ALO!  THEN  ALO!  = 
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AMP ! (C%) 


IF  AMP ! (C%)  >  AHI !  THEN  AHI !  =  AMP!(C%) 

C%  =  C%  +  1 

LOOP 

IF  K%  =  2  THEN 

FOR  1%  =  0  TO  C%  -  1 
J%  =  1%  +  1 

IF  J%  =  128  THEN  J%  =  127 
AMOP! (R%, 1%)  =  AMP! (J%) 

PROP! { R% , 1% )  =  PH! <J%) 

NEXT  1% 

K%  =  0 
ELSE 

FOR  1%  =  0  TO  C%  -  1 

AMOP ! { R% , 1% )  =  AMP ! ( I % ) 

PHOP! (R%, 1%)  =  PH! (1%) 

NEXT  1% 

END  IF 

CLOSE  #1 
NEXT  R% 

A$  =  "AMP" 

B$  =  "PHS" 

AHI!  =  AHI!  -  ALO! 

CALL  RANDOMIO { AMOP ! { ) , NN% , MM% , ALO ! , AHI ! , A$ ) 

CALL  ANGLEFIX ( PHOP ! ( ) , NN% , MM% ) 

FOR  R%  =  0  TO  MM% 

FOR  C%  =  0  TO  NN% 

IF  PHOP ! (R% , C% )  <  BLO!  THEN  BLO!  =  PHOP!(R%,C%) 

IF  PHOP ! (R% , C%)  >  BHI!  THEN  BHI !  = 

PHOP! (R% , C%) 

NEXT  C% 

NEXT  R% 

BHI!  =  BHI!  -  BLO! 

CALL  RANDOMIO ( PHOP ! ( ) , NN% , MM% , BLO ! , BHI ! , B$ ) 

PRINT  "***********  ALL  DONE  ************** 

END 

SUB  RANDOMIO  ( AMOP ! ( 2 ) , NN% , MM% , ZLO ! , ZHI ! , AS ) 

REM  THIS  SUBROUTINE  PRODUCES  RANDOM  ACCESS  FILES  OF  BOTH 
REM  THE  AMPLITUDE  DATA  AND  PHASE  DATA 
LOCAL  C% , R% , J% 

SHARED  P$ 

OPEN  "E: "+P$+A$+" .GRD"  AS  #2  LEN  =  4 
FIELD  #2,  4  AS  VERIFYS 
LSET  VERIFYS  =  "DSPM" :  PUT  #2 
FIELD  #2,  2  AS  BYTE1S ,  2  AS  BYTE2$ 

LSET  BYTE1S  =  MKI${100) 

LSET  BYTE2S  =  MKI$(64)  :  PUT  #2 

FIELD  #2,  4  AS  SP$ 

LSET  SP$  =  MKMS$ ( 0 ) :  PUT  #2 
LSET  SP$  =  MKMS$ (74):  PUT  #2 
LSET  SP$  =  MKMS$ ( 0 ) :  PUT  #2 


LSET  SP$  =  MKMS$ (51.1875) :  PUT  #2 
LSET  SP$  =  MKMS$ ( 0 ) :  PUT  #2 
LSET  SP$  =  MKMS$ ( ZHI ! ) :  PUT  #2 
FOR  R%  =  0  TO  MM% 

FOR  C%  =  0  TO  99 
J%  =  C%  +  14 

LSET  SP$  =  MKMS$(AMOP! (R%, J%)-  ZLO!)  :  PUT  #2 
NEXT  C% 

NEXT  R% 

CLOSE  #2 
END  SUB 

SUB  ANGLEFIX ( PHOP ! (2) , NN%,MM%) 

REM  THIS  SUBROUTINE  REMOVES  THE  EFFECTS  OF  THE  RECEIVER 
MODULO 

REM  360  DEGREE  EFFECT  ON  THE  DATA 
LOCAL  C%,  K%  ,  B ! 

REM  **************  WORKING  ALONG  THE  ROWS  ************** 
REM  THIS  PART  STARTS  AT  0 , 0  AND  MOVES  ALONG  THE  ROWS  TO 
MM , NN  FOR  R%  =  0  TO  MM% 

C%  =  0 

WHILE  C%  <  64 

IF  PHOP! (R%,C%+1)  -  PHOP ! (R% , C%)  >225  THEN 
FOR  K%  =  0  TO  C% 

PHOP ! (R% , K%)  =  PHOP ! (R% , K%)  +  360 
NEXT  K% 

END  IF 
INCR  C% 

WEND 

WHILE  C%  <  NN% 

IF  PHOP ! ( R% , C% )  -  PHOP! (R%,C%+1)  >225  THEN 
FOR  K%  =  C%+1  TO  NN% 

PHOP ! (R% , K% )  =  PHOP ! ( R% , K% )  +  360 
NEXT  K% 

END  IF 
INCR  C% 

WEND 
NEXT  R% 

REM*  ************************************************ 

REM* **********  WORKING  UP  THE  COLUMNS  ************** 

REM  THIS  PART  STARTS  AT  0,0  AND  MOVES  UP  THE  COLUMNS  TO 
MM , NN  FOR  C%  =  0  TO  NN% 

R%  =  32 
WHILE  R%  >  0 

WHILE  PHOP ! (R% , C%)  -  PHOP ! (R%-1 , C%) >200 

PHOP! (R%-1,C%)  =  PHOP! (R%-1,C%)  +  360 

WEND 

WHILE  PHOP! (R%-1 ,C%)  -  PHOP ! (R% , C%) >200 

PHOP! (R%-1,C%)  =  PHOP! (R%-1,C%)  -  360 

WEND 

INCR  R% , -1 


NEXT  C% 

FOR  C%  =  0  TO  NN% 

R%  =  32 

WHILE  R%  <  MM% 

WHILE  PHOP I { R% , C% )  -  PHOP! (R%+1,C%) >200 

PHOP! (R%+1,C%)  =  PHOP! (R%+1,C%)  +  360 

WEND 

WHILE  PHOP! (R%+1,C%)  -  PHOP !{ R% , C% }> 200 

PHOP! <R%+1,C%)  =  PHOP! (R%+1,C%)  -  360 

WEND 
INCR  R% 

WEND 
NEXT  C% 

REM  **************  WORKING  UP  THE  COLUMNS  ************* 
REM  THIS  PART  STARTS  AT  0 , 0  AND  MOVES  UP  THE  COLUMNS  TO 
MM , NN  FOR  R%  =  32  TO  MM% 

C%  =  27 
WHILE  C%  >  0 

WHILE  PHOP ! (R% , C%)  -  PHOP ! (R% , C%-1 )> 200 

PHOP! (R%,C%-1)  =  PHOP! (R%,C%-1)  +  360 

WEND 

WHILE  PHOP ! (R%,C%-1)  -  PHOP !( R% , C% )> 200 

PHOP! (R%,C%-1)  =  PHOP! (R%,C%-1)  -  360 

WEND 

INCR  C%,~1 

WEND 
NEXT  R% 

FOR  R%  =  32  TO  MM% 

C%  =  27 
WHILE  C%  <  64 

WHILE  PHOP ! { R% , C% )  -  PHOP !( R% , C%+1 )> 200 

PHOP! (R%,C%+1)  =  PHOP ! (R%,C%+1)  +  360 

WEND 

WHILE  PHOP ! (R%,C%+1)  -  PHOP !( R% , C% )> 200 

PHOP! (R%,C%+1)  =  PHOP! (R%,C%+1)  -  360 

WEND 
INCR  C% 

WEND 
NEXT  R% 

REM  ************************************************** 
FOR  R%  =  32  TO  MM% 

C%  =  90 
WHILE  C%  >  63 

WHILE  PHOP ! { R% , C% )  -  PHOP !( R% , C%-1 )> 200 

PHOP! (R%,C%-1)  =  PHOP ! (R%,C%-1)  +  360 

WEND 

WHILE  PHOP ! (R%,C%-1)  -  PHOP !( R% , C% )> 200 

PHOP! (R%,C%-1)  =  PHOP! (R%,C%-1)  -  360 


WEND 


INCR  C%,-1 


WEND 
NEXT  R% 

FOR  R%  =  32  TO  MM% 

C%  =  90 

WHILE  C%  <  NN% 

WHILE  PHOP ! (R% , C%)  -  PHOP !( R% , C%+1 )> 200 

PHOP ! ( R% , C%+1 )  =  PHOP i ( R% , C%+1 )  +  360 

WEND 

WHILE  PHOP! (R%,C%+1)  -  PHOP ! (R% , C% )> 200 

PHOP! (R%,C%+1)  =  PHOP!  (R%,C%+1)  -  360 

WEND 
INCR  C% 

WEND 
NEXT  R% 

REM  THIS  PART  STARTS  AT  0 , 0  AND  MOVES  UP  THE  COLUMNS  TO 
MM , NN  FOR  R%  =  0  TO  32 
C%  =  64 
WHILE  C%  >  0 

WHILE  PHOP ! ( R% , C% )  -  PHOP ! (R% , C%-1 ) >200 

PHOP! (R%,C%-1)  =  PHOP! (R%,C%-1)  +  360 

WEND 

WHILE  PHOP! (R%,C%-1)  -  PHOP ! (R% , C%) >200 

PHOP! (R%, C%-1)  =  PHOP! (R%,C%-1)  -  360 

WEND 

INCR  C%,-1 

WEND 
NEXT  R% 

FOR  R%  *  0  TO  32 
C%  =  64 

WHILE  C%  <  NN% 

WHILE  PHOP ! { R% , C% )  -  PHOP ! (R% , C%+1 )> 200 

PHOP! (R%,C%+1)  =  PHOP ! (R%,C%+1)  +  360 

WEND 

WHILE  PHOP! (R%,C%+1)  -  PHOP ! (R% , C%) > 200 

PHOP! (R%,C%+1)  =  PHOP ! (R%, C%+1)  -  360 

WEND 
INCR  C% 

WEND 
NEXT  R% 

END  SUB 
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Appendix  E.  Fast  Fourier  Transform  Program 


PRINT  "**************  FFTAZEL . BAS  ***************  » 

PRINT  "THIS  PROGRAM  IS  DESIGNED  TO  TAKE  DATA  FROM  A  2-D  ” 
PRINT  "MATRIX  ARRAY  REPRESENTING  THE  NEAR  FIELD  IN  " 

PRINT  "AMPLITUDE  (DB)  AND  PHASE  (DEG)  OF  AN  ANTENNA  AND  " 
PRINT  "TURN  IT  INTO  THE  FAR  FIELD  IN  AMPLITUDE  AND  PHASE  " 
PRINT  "OF  THE  ANTENNA  IN  THE  FAR  FIELD. 

PRINT  "THE  OUTPUT  IS  CONVERTED  TO  AZIMUTH  AND  ELEVATION" 
PRINT  "VALUES  FOR  EASY  COMPARISON  WITH  MEASURED  VALUES." 
PRINT  "PLEASE  ENSURE  THAT  THE  DATA  DISKETTE  IS  IN  DRIVE  A:" 
PRINT  "ENSURE  DATA  FILES  ARE  LABELLED  0-360FDT1-64 . BAS . " 
PRINT  "RESULTS  FOR  PLOTTING  WILL  BE  IN  E : 0-360FFA . GRD" 

PRINT  "HOW  MANY  DEGREES  ARE  THE  PLATES  ABOVE  THE  DISH 
INPUT  P% 

INPUT  "HOW  MANY  ROWS  IN  THE  MATRIX  TO  BE  PROCESSED  " , M% 
MM%=M%-1 

mu%  =  cint ( log (m% ) /log ( 2 ) ) 

INPUT  "HOW  MANY  COLUMNS  IN  THE  MATRIX  TO  BE  PROCESSED  " , N% 
NN%=N%- 1 

nu%  =  cint (log (n%) /log (2) ) 
pi#  =  3.141592653589793 

DIM  XR ! ( 0 : MM% ,  0 : NN% ) , XI ! ( 0 : MM% ,  0 : NN% ) 

DIM  TXR ! ( 0 : MM% ,  0 : NN% ) , TXI ! ( 0 : MM% ,  0 : NN% ) 

DIM  AMOP! (0:MM%,  0 : NN% ) , PHOP 1(0: MM% ,  0:NN%) 

DIM  XXR ! (0:NN%) ,  XXI!(0:NN%) 

DIM  XXRR! (0:MM%) ,  XXIR!(0:MM%) 

DIM  AMP ! ( 0 : NN% ) ,  PH ! ( 0 : NN% ) ,  PHASE ! ( 0 : NN% ) 

PRINT  "***************  WORKING  *****************" 

FOR  I%=0  TO  MM% 

IF  1%  <  9  THEN 

B$  =  RIGHT$(STR$(I%+1) ,1) 

ELSE 

B$  =  RIGHTS (STR$ <I%+1) , 2) 

END  IF 

IF  P%  <  10  THEN 

P$  =  RIGHTS (STR$ (P%) , 1) 

ELSEIF  P%  >  9  AND  P%  <  100  THEN 
P$  =  RIGHTS (STR$ (P%) , 2) 

ELSE 

P$  =  RIGHTS (STR$ (P%) , 3) 

END  IF 

c  $  =  "A: "  +  P$  +  "FDT"+B$  +  " .BAS" 
open  c$  for  input  as  #1 
FOR  JL%=0  TO  NN% 

input  #1,  AMP!  ( j 1 % )  ,PH!  (  j 1 % ) 

PHASE!  ( j 1% )  =  2*pi#*PH!  (jl%)/360 

XR!  (1%, JL%)  =10'  (AMP!  (JL%) / 1 0 )  *COS (PHASE!  (JL%)  ) 
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XI!  (1%, JL%)  =10  * (AMP!  { JL% ) / 10 ) *  SIN (PHASE!  (JL%)  ) 
NEXT  JL% 

CLOSE  #1 

NEXT  1% 

FOR  R%  =  0  TO  MM% 

FOR  C%  =  0  TO  NN% 

::XR!  (C%)  =  XR! (R%,C%) 

XXI ! (C%)  =  XI ! ( R% , C% ) 

NEXT  C% 

CALL  FFT ( XXR ! ( ) , XXI ! ( ) , N% , nu% ) 

FOR  C%  =  0  TO  NN% 

XR! (R% , C%)  =  XXR! (C%) 

XI! (R%,C%)  =  XXI! (C%) 

NEXT  C% 

NEXT  R% 

REM 

******************************************************** 
REM  *  *  *  *THIS  PORTION  OF  THE  PROGRAM  DOES  A  CUT  AND***** 
REM  **********  DIAGONAL  SWAP  OF  THE  DATA  VALUES.******* 
REM*  ***********************************  *  *************** 
FOR  C%  =  0  TO  NN% 

FOR  R%  =  0  TO  MM% 

XXRR ! (R%)  =  XR! (R%, C%) 

XXIR ! (R%)  =  XI ! (R%,C%) 

NEXT  R% 

CALL  FFT ( XXRR ! ( ) , XXIR ! ( ) , M% , mu% ) 

FOR  R%=0  TO  MM% 

XR! (R%,C%)  =  XXRR! (R%) 

XI! ( R% , C% )  =  XXIR! (R%) 

NEXT  R% 

NEXT  C% 

PRINT  "***************  STILL  WORKING  ***************" 
FOR  R%  =  0  TO  MM% 

FOR  C%  =  0  TO  NN% 

J%  =  C%  +  N%/2 

IF  J%  >  NN%  THEN  J%  =  J%  -  N% 

TXR ! ( R% , C% )  =  XR ! ( R% , J% ) 

TXI ! ( R% , C% )  =  XI! (R%, J%) 

NEXT  C% 

NEXT  R% 

FOR  C%  =  0  TO  NN% 

FOR  R%  =  0  TO  MM% 

K%  =  R%  +  M%/2 

IF  K%  >  MM%  THEN  K%  =  K%  -  M% 

XR ! (R% , C%)  =  TXR ! ( K% , C% ) 

XI ! (R% , C%)  =  TXI ! ( K% , C% ) 

NEXT  R% 

NEXT  C% 

PRINT  "**************  STILL  AT  IT  ****************" 
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M-UM  943 
UNCLASSIFIED 


HOOIFICATION  OF  PARABOLIC  DISH  ANTENNA  PATTERN  USING 
TWO  SVNAETRICALLV  PL.  .  (U>  AIR  FORCE  INST  OF  TECH 
HRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGI. .  G  C  THORPE 
DEC  G7  AFIT/GE/ENQ/87D-67  F/6  9/1 


2/2 


,"*<>*~*,‘,‘ *'*'*'*  **♦  *'♦  *'*  ***•*'*  **>-**»■♦  *■*  *•'  *•*  ***  ***-*** 


1.0  u  •-  la  1 
==  “  I-.  *2.2 


M_  L“  la  i 


1.8  i 


L25  I U  in  i. 


ZLO !  =  9 . 9E+30 
ZHI!  =  -ZLO ! 

FOR  R%  =  0  TO  MM% 

FOR  C%  =  0  TO  NN% 

AMOP !  (R% , C%)  =  10*LOG ( SQR (XR !  (R% , C%) ^  2+XI !  (R% , C%) *  2 ) ) 
IF  AMOP ! (R% , C%)  <  ZLO!  THEN  ZLO!  =  AMOP!(R%,C%) 

IF  AMOP ! (R% , C%)  >  ZHI!  THEN  ZHI!  =  AMOP!(R%,C%) 

IF  XR! (R% , C%)  =  0  THEN 
PHOP ! (R% , C%)  =  90 
ELSE 

PHOP! (R%  C%)  = 

( ATN ( XI !  (R% , C%) /XR !  <R%,C%) ) *  360 ) / ( 2*PI# ) 

END  IF 
NEXT  C% 

NEXT  R% 

CALL  DATAFILE (AMOP ! ( ) , NN% , MM% , ZLO ! , ZHI ! ) 

PRINT  "***********  ALL  DONE  *************" 

END 

SUB  FFT(XXR! (1) , XXI! (1) ,N%,NU%) 

SHARED  PI# 

LOCAL  J% , C% , KK% , K% , B , 1% , L%, K2% , TRE ! , TIM ! , XRT ! , XIT ! 

FOR  I%=0  TO  N%-1 
J%=0 
KK%=I% 

K%=1 

TEST: 

B*KK%/2~ (NU%-K%) 

IF  B  >■  1  THEN 

J%  *  J%  +  2~(K%  -  1) 

KK%=KK%  -  2~ (NU%  -  K%) 

END  IF 
K%=K%+1 

IF  K%  >  NU%  THEN 
GOTO  CHK 
ELSE 

GOTO  TEST 
END  IF 

CHK: 

IF  J%  <  1%  THEN 
GOTO  JMP 


ELSE 


TRE!=XXR! (J%) 
TIM ! =XXI ! (J%) 
XXR ! (J%)=XXR! (1%) 
XXI! (J%)=XXI! (1%) 
XXR! ( 1% ) =TRE ! 

XXI! (1%) =TIM! 
END  IF 


JMP: 


NEXT  1% 


FOR  L%=0  TO  NU%-1 
C%=1 

K2%=2~L% 

P#=PI#/K2% 

FOR  K%=0  TO  N%-1 

XRT ! =XXR ! (K%+K2%) *COS (K%*P#) +XXI ! (K%+K2%) *SIN(K%*P#) 

XIT ! =XXI ! (K%+K2%) *COS(K%*P#) -XXR! (K%+K2%) *SIN(K%*P#) 

XXR! ( K%+K2% ) =XXR ! (K%)-XRT! 

XXI! (K%+K2%) =XXI ! (K%)-XIT! 

XXR! (K%) =XXR ! (K%)+XRT! 

XXI! (K%)=XXI! (K%) +XIT! 

IF  C%=K2%  THEN 
K%=K%+K2% 

C%=0 
END  IF 
C%=C%+1 
NEXT  K% 

NEXT  L% 

END  SUB 

***********************  DATA  FILE  ******************* 
SUB  DATAFILE  ( AMOP ! ( 2 ) , NN% , MM% , ZLO ! , ZHI ! ) 

LOCAL  C%,R% 

SHARED  P$,PI# 

AS  *  "A: "+P$+"FFAA. DAT” 

B$  *  “A: "+P$+"FFKS • DAT" 

OPEN  A$  FOR  OUTPUT  AS  #1 
OPEN  B$  FOR  OUTPUT  AS  #2 
1%  =  -9732 
FOR  R%  =  18  TO  45 
J%  =  -9995 

I!  =  IV10000  :REM  Ky  VALUE 

FOR  C%  =  42  TO  85 

J!  -  J%/10000  :REM  Kx  VALUE 

ELL!  =  ATN ( I ! / SQR ( 1-1 ! ~  2 ) ) 

AZARG!  =  J ! /COS ( ELL ! ) 

IF  1 -AZARG ! *  2  <  0  THEN  GOTO  XMP 
AZZ !  =  ATN (AZARG ! /SQR ( 1- AZARG ! ~  2 ) ) 

AZH!  =  (AZZ ! /PI#) *180 
ELV!  =  (ELL! /PI#) *180 

IF  AMOP ! (R% , C%)  <  30  THEN  AMOP!(R%,C%)  =  30 
Q$  =  LEFTS (STR$ (J! ) , 9) 

R$  =  "  "+LEFT$ ( STR$ (I ! ) , 9) +" 

S$  =  LEFTS (STR$( AMOP ! (R%fC%) ) ,8) 

PRINT  #2,  Q$ , R$ , S$ 

AMOP ! (R% , C%)  =  AMOP ! (R% , C%) *COS (AZZ ! ) *COS (ELL ! ) 

Q$  =  LEFTS (STR$ (AZM! ) , 9) 

R$  ■  "  "+LEFTS (STR$ (ELV! ) , 9) +" 

S$  =  LEFTS (STR$ (AMOP ! <R%,C%) ) ,8) 


:REM  Ky  VALUE 


:REM  Kx  VALUE 


1 

1 


-  V.V  “ 


/  »  »  •  •  « 


“  ■  •  ■  ,  ■  "Jfc  "  » 
*  -  •  v" 


^  i 

v, 


PRINT  #1,  Q$ ,R$ , S$ 
XMP : 

INCR  J% , 464 
NEXT  C% 

INCR  1%,  721 
NEXT  R% 

CLOSE  #1 
CLOSE  #2 
END  SUB 
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PRINT  "THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  DISTANCES" 

PRINT  "TO  THE  REFERENCE  PLANE  AND  THE  X  CO-ORDINATES" 

PRINT  "OF  THE  POINTS  THAT  REFLECT  OFF  THE  MOVEABLE  PLATE" 
PRINT  "THE  PLATE  IS  MOVED  OUT  1/16  WAVELENGTH  EACH  TIME" 
PRINT  "FROM  POSITION  0  TO  POSITION  8." 

PRINT  "WHERE  0  MEANS  AT  THE  LOWEST  POSITION  AND  8  IMPLIES" 
PRINT  "THAT  THE  PLATE  IS  1/2  WAVELENGTH  ABOVE  THE  SURFACE" 
INPUT  "WHAT  POSITION  IS  THE  PLATE  IN?  " ; L% 

REM  FROM  THIS  VALUE  OF  L%  THE  HEIGHT  IN  INCHES  IS  FOUND 

L ! =L%* • 230625+ . 2 

L$  =  RIGHTS ( STR$ ( L% ) , 1 ) 

REM  TPI !  =  2  PI 
TPI!=2*3. 141592654 

REM  B!  IS  THE  ANGLE  THE  TANGENT  MAKES  WITH  THE  X-AXIS 
REM  THIS  ANGLE  HAS  A  VALUE  OF  14.6  DEGREES 
B!  =  0.25481876 

REM  TH!  IS  THE  ANGLE  MEASURED  CLOCKWISE  FROM  THE  X-AXIS 
REM  TO  THE  RAY  THAT  INTERSECTS  WITH  THE  TANGENT  AND  THE 
REM  PARABOLA  ; ITS  VALUE  IS  60.8  DEGREES 
TH!=1. 0611588 

REM  FIRST  STEP  IS  TO  FIND  THE  CENTER  CO-ORDINATES  OF  THE 
PLATE 

REM  8.986  IS  THE  X  CO-ORDINATE  ON  THE  PARABOLA  SURFACE 

REM  16.08  IS  THE  Y  CO-ORDINATE  ON  THE  PARABOLA  SURFACE 

XC ! =8 . 986-L ! *SIN (B ! ) 

YC ! =16 . 08-L ! *COS (B! ) 

REM  PLATE  CENTER  X  CO-ORD  =  XC ! 

REM  PLATE  CENTER  Y  CO-ORD  =  YC ! 

REM  OPEN  FILES  TO  ACCEPT  AMPLITUDE  AND  PHASE  VALUES 
CALCULATED 

A$  =  "A:PLAMP"+L$+" .DAT" 

B$  =  "A:PLPHS"+L$+" .DAT" 

OPEN  A$  FOR  OUTPUT  AS  #1 
OPEN  B$  FOR  OUTPUT  AS  #2 

REM  NOW  FIND  THE  LEFT  EDGE  CO-ORDINATES  OF  THE  PLATE 
REM  PRC!  IS  THE  CO-ORDINATE  ON  THE  SURFACE  OF  THE  PLATE 
FOR  PRC! =-2.5  TO  2 . 5  STEP  .1 

REM  XL!  IS  THE  X  CO-ORDINATE  ON  THE  SURFACE  OF  THE  PLATE 
XL ! =XC ! +PRC ! *COS ( B ! ) 

REM  YL!  IS  THE  Y  CO-ORDINATE  ON  THE  SURFACE  OF  THE  PLATE 
YL ! =YC ! -PRC ! *  SIN { B ! ) 

IF  PRC! =-2.5  THEN  XPL!=XL! 

REM  X  CO-ORD.  OF  LEFT  MOST  POINT  ON  PLATE 
REM  FIND  ANGLE  FROM  FOCUS  TO  THE  EDGE  OF  PLATE 
PHI#=ATN{YL! /XL! ) 

IF  PRC!  =  -2.5  THEN  PHIL#=PHI# 


REM  PHIL*  IS  THE  ANGLE  TO  THE  LEFT  EDGE  OF  THE  PLATE 
IF  PRC  !  >  2.4  THEN  PHIR#=PHI# 

REM  PHIR*  IS  THE  ANGLE  TO  THE  RIGHT  EDGE  OF  THE  PLATE 
REM  FIND  CW  COMP.  FROM  THE  VERT.  OF  REFLECTION  ANGLE 
VC#»1 . 5707 963 27-PHI #-2*B ! 

REM  FIND  THE  DISTANCE  FROM  THE  FOCUS  TO  THIS  POINT 
DTLE!«SQR(XL! ~2+YL! ~2) 

REM  FIND  DISTANCE  TO  THE  REFERENCE  PLANE 
DTRP ! = ( 12 . 25+YL ! ) /COS (VC# ) 

REM  FIND  TOTAL  DISTANCE  FROM  FOCUS  TO  REFERENCE  PLANE 
TOTD ! =DTRP ! +DTLE ! 

REM  FIND  X  CO-ORDINATE  OF  RAY  IMPACT  IN  REF.  PLANE 
XCRP ! * ( (TAN (VC#) ) * ( YL ! +12 . 25 ) )+XL! 

IF  PRC!  =  -2.5  THEN  LXCRP ! =XCRP ! 

REM  LXCRP!  IS  THE  LEFT  MOST  CO-ORD  OF  REFLECTION  ZONE 
IF  PRC!  >  2.4  THEN  RXCRP ! =XCRP ! 

REM  RXCRP!  IS  THE  RIGHT  MOST  CO-ORD  OF  REFLECTION  ZONE 
NEXT  PRC! 

REM  LTPRL!  IS  THE  DISTANCE  FROM  THE  FOCUS  TO  THE  PARABOLA 

REM  OF  THE  RAY  THAT  JUST  MISSES  THE  LEFT  EDGE  OF  THE  PLATE 

LTPRL ! =34. 5/ ( 1-SIN ( (2*3. 141592654) -PHIL#) ) 

XCPRL ! =LTPRL ! *COS ( PHIL# ) 

REM  XCPRL!  IS  CO-ORD  OF  POINT  ON  PARABOLA 
LTPRR ! =34. 5/ (1-SIN ( (2* 3. 141592654) -PHIR#) ) 

REM  LTPRR !  IS  THE  DISTANCE  FROM  THE  FOCUS  TO  THE  PARABOLA 

REM  OF  THE  RAY  THAT  JUST  MISSES  THE  RIGHT  EDGE  OF  THE  PLATE 

XCPRR ! =LTPRR ! *COS ( PHIR# ) 

REM  TO  THE  RIGHT  OF  XCPRR!  DIRECT  RAYS  CAN  REACH  THE 
REF.  PLANE 

REM  FIND  THE  AMPLITUDE  OF  THE  REFLECTED  ENERGY 
REAMP ! =  (XCPRR ! -XCPRL ! ) / (RXCRP ! -LXCRP ! ) 

REM  REPEAT  ABOVE  LOOP  WITH  KNOWN  VALUES 
FOR  PRC! =-2.5  TO  2.5  STEP  .1 
XL! =XC ! +PRC ! *COS (B! ) 

YL ! =YC ! -PRC ! *SIN (B ! ) 

REM  FIND  ANGLE  FROM  FOCUS  TO  THE  EDGE  OF  PLATE 
PHI#=ATN ( YL ! /XL ! ) 

REM  FIND  CW  COMP,  FROM  THE  VERT.  OF  REFLECTION  ANGLE 
VC#=1 . 570796327-PHI#-2*B ! 

REM  FIND  THE  DISTANCE  FROM  THE  FOCUS  TO  THIS  POINT 
DTLE ! =SQR (XL ! ~2+YL! ~2) 

REM  FIND  DISTANCE  TO  THE  REFERENCE  PLANE 
DTRP ! = ( 12 . 25+YL ! ) /COS ( VC# ) 

REM  FIND  TOTAL  DISTANCE  FROM  FOCUS  TO  REFERENCE  PLANE 
TOTD ! =DTRP ! +DTLE ! 

TOTD ! =46 . 7  5-TOTD ! 

PHSR ! = ( (TOTD 1/3.69) *TPI ! ) 

REM  FIND  X  CO-ORDINATE  OF  RAY  IMPACT  IN  REF.  PLANE 
XCRP ! = ( (TAN (VC#) ) * (YL1+12 . 25) ) +XL! 

XPOS ! =64+ (XCRP 1/39.685) *64 


IF  PHSR!  >  1.5707963  THEN  PHSR!  =  PHSR!  -  1.5707963 
IF  XCRP I <  XPL !  THEN 

RFPXA ! *0.5* (COS (PHSR!) )  +  1 
RFPYA ! =0 . 5 * SIN { PHSR !) 

RFPA!  =SQR(RFPXA!~2+RFPYA!~2) 

PHSE!  =ATN (RFPYA! /RFPXA ! ) 

ELSEIF  XCRP ! >XCPRR !  THEN 

RFPXA! =0.5* (COS (PHSR! ) )  +  1 
RFPYA! =0.5* SIN (PHSR! ) 

RFPA!  =SQR  (RFPXA !  ""2+RFPYA !  * 2 ) 

PHSE!  =ATN (RFPYA! /RFPXA ! ) 

ELSE 

RFPXA !=0. 5*COS (PHSR! ) 

RFPYA ! =0.5* SIN ( PHSR ! ) 

RFPA!  =SQR(RFPXA! ~2+RFPYA! *2) 

PHSE!  =ABS (ATN (RFPYA! /RFPXA! ) ) 

END  IF 

PHSE!  =( (PHSE!/(TPI! ) ) *360) 

WRITE  #1,  XPOS ! , RFPA ! 

WRITE  #2,  XPOS!, PHSE! 

NEXT  PRC! 

CLOSE  #1 
CLOSE  #2 
END 
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